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Preface 

It  is  the  purpose  of  this  thesis  to  describe  several  Methods  of 
obtaining  neutron  flux-spectra  information  from  foil  data*  These 
methods,  although  presented  In  considerable  detail,  are  not  jet  in 
optimized  fora*  It  is  hoped  that  this  can  be  done  in  the  near  future* 
lb  Is  veil  known  that  threshold  and  fission  type  foils  are  most 
advantageously  used  to  generate  integral  flux  data.  Since  radiation 
effects  investigators  usually  desire  dose  parameters  rather  than  flux 
parameters,  a  method  of  calculating  the  dose  from  integral  flux  data 
is  presented  in  Appendix  A* 

The  author  wishes  to  thank  Dr.  V.  J.  Price  and  Lt.  Col.  B.  Z. 
Robertson  for  their  helpful  advice  and  consideration*  The  author  is 
greatly  indebted  to  Lt.W.  R*  Burrus,  Materials  Laboratory,  for  his 
many  pertinent  suggestions  and  comments.  Dr.  M.  A.  Edvards,  Lt.  T.  A* 
Brown,  and  their  associates  of  the  Computations  Brandi,  Aeronautical 
Research  Laboratory,  suggested  some  of  the  mathematical  techniques 
and  handled  the  programming  of  the  digital  computer.  Mr.  S.  A* 
Szavlevicz,  Materials  Laboratory,  gratiously  made  the  necessary  finan¬ 
cial  arrangements  for  the  computer  program.  Without  the  ideas  and 
assistance  of  any  of  the  above  mentioned  people,  this  thesis  vould 
not  have  been  possible. 

Paul  M.  Utbe,  Jr. 
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Abstract 

Several  analytical  methods  are  described  for  obtaining  neutron 
flux-spectra  knowledge  from,  foil  data*  The  response  functions  used 
ere  the  nuclear  reaction  characteristics  of  threshold*  fast  fission* 
and  resonance  type  foils* 

Two  previously  used  methods  and  two  nei'  methods  are  discussed* 

The  previously  used  methods  are  the  Trice  Method  and  tbs  Cadmium 
Difference  Method*  The  Trice  Method*  applicable  in  tbe  fast  neutron 
energy  region,  uses  a  *step  function"  cross  section  representation* 

The  step  function  energy  threshold  is  determined  by  assuming  a  differ* 
ential  flux  function*  Tbe  Cadmium  Difference  Method*  applicable  in 
the  l/E  energy  region,  uses  the  resonance  activation  integrals  of 
resonance  type  foils  as  cross  section  parameters*  Although  consid¬ 
eration  is  given  to  these  methods*  major  interest  is  focused  on  two 
new  methods*  These  methods  are  called  the  Polygonal  Method  and  the 
Polynomial  Method* 

The  Polygonal  Method  uses  the  cross  sections  in  their  true  fora 
and  repreeents  the  flux-energy  spectrum  by  a  combination  of  linear 
functions  of  energy*  In  performing  the  necessary  calculations*  diffi¬ 
culty  in  selecting  certain  energy  parameters  is  experienced* 

The  Polynomial  Method  uses  the  cross  sections  is  their  true 
form  and  represents  the  flux-energy  spectrum  by  s  weighting  function 

Ti 
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timer  a  polynomial*  The  degree  of  the  polynomial  la  one  less  than 
the  number  of  different  foils  used*  A  digital  computer  is  conveniently 
used  for  the  required  numerical  calculations* 

The  result  of  several  test  cases  using  the  Polynomial  Method 
indicates  that  considerable  improvement  over  the  Tfcice  Method  is 


achieved, 


ATTAINMENT  (F  NEUTRON  FIUT-SIECTRA 
FROM 

FOIL  ACTIVATIONS 


I.  Introduction 


Tbs  purpoas  of  this  thesis  is  to  describe  several  methods  for 
obtaining  neutron-flux  spectra*  in  high  intensity  reactors.  Know¬ 
ledge  of  the  neutron  flux  environment  is  required  for  the  effective 
interpretation  of  the  results  obtained  from  radiation  effects  exper¬ 
iments.  The  analysis  methods,  described  are  designated  ast 

A.  The  Trice  Method 

B.  The  Cadmium  Difference  Method 

C.  The  Polygonal  Method 

D.  The  Polynomial  Method 

These  methods  are  presented  In  Sections  3*  5>  and  6  respectively. 

The  analysis  methods  are  based  upon  the  nuclear  reaction  charac¬ 
teristics  of  three  types  of  foils;  (l)  threshold,  (2)  fast  fission, 
and  (3)  resonance.  The  word  ■foil"  does  not  imply  a  physical  charac¬ 
teristic;  i.e.,  a  foil  may  be  in  powder,  pellet,  liquid  or  other 

*  In  the  remainder  of  this  report,  "spectra"  signifies  "energy  spec¬ 
tra."  ’ 
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physical  form.  The  word  ■radloactivant'*,  though  not  commonly  used, 
would  be  more  appropriate. 

Major  Interest  is  focused  on  the  threshold  and  fast-fission  types 
of  foils.  For  tbs  sake  of  generalty,  some  brief  comments  regarding 
resonance  type  foils  are  included.  It  is  be  lie-red,  however,  that  the 
calculational  methods  derived  are  sufficiently  general  so  as  to  be 
applicable  to  any  detector  system. 

Of  the  many  techniques  that  have  been  used  to  determine  reactor 
neutron  spectra,  the  foil  technique  is  the  most  often  used.  This  is 
so,  not  because  the  foil  gives  the  most  precise  information,  but 
because  it  has  the  greatest  compatibility  with  the  nuclear  and  phys¬ 
ical  reactor  characteristics.  Three  specific  reasons  for  this  general 
acceptance  are  as  follows  *■  (l)  only  a  few  milligrams  of  foil  are 

needed  in  most  cases;  (2)  the  foil,  with  suitable  cover,  is  the  com¬ 
plete  detector;  and  (3)  all  analysis  is  done  in  a  convenient  nuclear 
counting  room  after  irradiation.  It  is  not  surprising  that  the 
radioactivant  foil  has  achieved  such  popularity  in  this  era  of  large 
and  complicated  experiments. 

On  the  otber  hand,  the  use  of  foils  does  introduce  several  com¬ 
plications.  Extreme  care  must  be  taken  in  the  processing  of  the 
foils  (both  before  and  after  irradiation).  If  this  is  not  done, 
misleading  resulte  may  be  obtained.  In  the  end,  the  proper  use  of 
foils  reduces  to  an  art  rather  than  a  science,  for  once  the  correct 


irradiation  and  counting  data  are  recorded,  spectral  infcntattcn  cm  be 
generated  from  the  original  data  as  better  analytical  methods  become 
available.  In  this  study,  analytical  methods  trill  be  emphasized.  As 
of  this  writing,  the  msthods  have  not  been  subjected  to  a  t>ood  test 
or  intercomparisonj  however,  it  is  hoped  that  this  can  be  done  in  the 


near  future 
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IX*  Characteristics  of  foils 


General 

Before  the  analytical  methods  ean  be  discussed,  it  is  necessary 
to  discuss  the  nuclear  characteristics  of  certain  types  of  foils* 
This  is  necessary  because  the  usefulness  of  the  analysis  methods  is 
directly  dependent  upon  the  neutron  reaction  characteristics  of 
certain  types  of  foils* 

In  general*  when  materials  are  Irradiated  in  a  neutron  flux, 
neutron  Induced  transmutation  of  their  nuclei  may  occur*  As  an 
example,  neutron-gamma  (n ,V).  neutron-proton  (n,p),  and  neutron- 
alpha  (n,0  reactions  occur  in  aluminum-27  when  this  material  is 
Irradiated  in  a  nuclear  reactor*  Bach  of  these  reactions  leads  to 
a  nev  and  different  isotope,  each  isotope  being  radioactive*  It  is 
the  radioactive  decay  of  the  isotopes  vhich  is  observed  in  the 
counting  room. 

The  reaction  probabilities  or  cress  sections  are  functions  of 
neutron  energy*  It  is  the  unusual  nature  of  the  reaction  cross 
sections  of  several  materials  vhich  is  responsible  for  their  use- 
f  .lro.sa  as  neutron  flux-spectra  detectors* 

The  number  of  particular  transmutations  occurring  in  a  unit 
interval  of  time,  per  unit  volume  of  material,  due  to  neutrons  having 
energies  in  the  energy  range  dX  about  X  1st 
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n(s)  d£  =  f{S)  tJ-(I)  N  dl 
^  «  -  <j>W  <£<*)  d£ 


Where:  R(E)  dB  «  the  number  of  transmutations  per  ce  per  see 
due  to  neutrons  In  the  energy  range  dS 
about  X. 

fi(c)  -  the  differential  neutron  flux  in  neutrons 
per  cn?  per  sec  per  energy  interval. 

-  the  probability  of  the  .ith  type  of  trans¬ 
mutation  occurring  per  unit  distance  of 
travel  of  a  neutron  with  energy  £. 

» 

N  -  the  number  of  pertinent  isotopic  nuclei 
per  cc. 

<Tc(e)  •=.  the  microscopic  cross  section  for  the  ith 
type  of  process  in  cm?.  As  mentioned 
previously,  the  microscopic  cross  section 
is  usually  expressed  in  units  of  barns. 

As  an  illustration  of  hov  a  particular  type  of  cross  section  can 
be  used  to  obtain  spectral  information,  the  following  example  is 
-  presented: 

Let  the  activation  cross  section  be  a  finite  constant  for  all 
neutron  energies  greater  than  and  less  than  £3  ^  zero  for  all 
other  energies.  A  graphical  plot  of  this  highly  idealized  cross 
section  is  given  in  Fig,  1, 
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Upon  Irradiation  in  an  unknown  flux,<j)  (E),  the  total  flux  between 
E^  and  Eg  can  b*  determined  as  the  integral  of  the  differential  flux. 
This  is  shown  bj  the  following* 


SOO  rEg 

d-(E)  <j>(E)  dE  =  Rfl^j  <j>(E)  dS 

Solving  for  the  total  flux  integrals 

2  $(2)  dE  =  ^ 

®1 

It  is  assumed  that  the  result  of  the  ith  type  of  transmutation  gives 
an  isotope  whose  quantity  at  any  time  can  be  determined,  that  the  foil 
did  not  perturb  the  flux,  and  that  the  flux  is  time  independent.  . 
Threshold  Type  Foils 

A  graphical  plot  of  the  Idealized  threshold  type  cross  section  Is 
shown  in  Fig.  ?. 

From  the  calculus,  a  result  similar  to  the  one  previously  obtained 
can  be  derived  as  shown  below** 


q(E)  dS 


(?) 


It  is  seen  that  this  type  of  foil  will  determine  the  total  neutron 
flux  above  Ej.  If  several  foils  are  used  which  have  differing  Ej*s, 
it  is  possible  to  determine,  in  theory,  the  total  flux  between  any  two 
thresholds.  In  this  manner,  a  histogram  type  of  spectra  representation 

I  « 
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*  A  =  -j  and  will  be  called  the  "saturated  activity* 
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can  be  generated  by  single  subtraction. 

Lhlortunately*  no  materials  exhibit  such  an  idealized  threshold. 
Figure  3  is  a  graphical  presentation*  with  curve  smoothing*  of  the 
cross  sections  of  some  of  the  cormonly  used  threshold  type  foils.  The 
fast  fission  foils  are  included  for  they  also  are  of  this  nature, 
indeed  even  more  so* 

The  variability  and  non-analytical  behavior  of  these  cross  sections 
are  responsible  for  the  difficulty  in  obtaining  good  spectral  infor¬ 
mation  from  foil  measurements. 

Resonance  Type  Foils 

If  the  foil  cross  section  shown  in  Fig.  1  were  of  such  a  width 
in  energy  that  (J>(E)  could  be  considered  a  constant  in  this  range* 
then  by  the  calculus i 


(V) 


A  =  ^°<|>(E)  q(B)dB  *  iq  (fCEj)  [e2  -  Ex] 

Where  <£( Ex)  =  <£(E2)  =  a  constant. 

By  this  method*  the  differential  flux  at  E^  can  be  obtained.  If 
several  foils  with  narrow  and  separated  reaction  intervals  (E^  to  E2) 
were  available,  the  differential  flux  at  the  pertinent  energies  could 
be  obtained.  Drawing  a  smooth  curve  through  the  plotted  points  would 
yield  a  reasonable  graph  of  the  differential  neutron  flux  (excluding 
fine  strr-ture).  If  materials  existed  with  these  cross  section  char¬ 
acteristics  and  if  their  reaction  intervals  were  sufficiently  separated 
30  that  the  total  energy  interval  would  be  1  ev  to  10  Mev,  then  the 
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problem  of  obtaining  a  neutron  flux-opectnaa  vould  bo  greatly  reduced-. 
Unfortunately,  this  Is  not  the  case. 

The  resonance  type  foils  that  have  been  used  in  the  past  have 
cross  sections  obeying,  in  general,  the  Briet-Wigner  relationship. 

Few  of  these  cross  sections  have  actually  been  measured  directly,  but 
seal-quantitative  values  can  be  calculated  by  using  the  Briet-Wigner 
formula  with  the  measured  energy  widths,  peak  energies,  and/or  maxi- 
aun  cross  sections  {jjil5j.  A  simpler  cross  section  parameter  is 
thP  empirically  determined  resonance  activation  integral}  however, 
reliance  on  resonance  integrals  often  leads  to  misleading  and  inaccu¬ 
rate  results  because  of  the  difficulty  in  taking  into  account  all 
the  necessary  corrections. 

For  a  more  complete  discussion  on  the  resonance  problem,  the 
reader  is  invited  to  read  Hughes'  Pile  Neutron  Research.  Suffice  it 
to  say  here  that  a  real  cross  section  problem  exists  in  the  resonance 
region,  that  is,  between  100  Kev  and  .02$  ev. 


in.  The  Trice  Method  and  the  Cadndm 


Difference  Method 


Introduction 

Present  evidence  indicates  that*  in  general,  any  reactor  neutron 

spectrum  above  the  thermal  energies  may  be  divided  into  two  parts. 

One  of  these  exists  at  energies  greater  than  1  Mev  and  has  the  general 

characteristics  of  a  fission  spectrum.  The  other,  a  result  of  neutrons 

being  "slowed-down"  hy  elastic  collisions,  occurs  at  energies  below 

1  Mev  and  has  an  approximate  1*^5  behavior.  Although  these  general 

postulates  are  crude,  they  have  been  used  with  notable  success. 

For  lack  of  a  more  appropriate  name,  the  fast  neutron  technique 

discussed  in  this  section  is  called  the  "Trice  Method".  Hr.  J.  B. 

Trice,  then  of  the  Oak  Ridge  National  Laboratory,  originally  instructed 

the  author  in  the  usefulness  pf  this  method. 

The  method  used  for  the  ife  spectral  region  is  commonly  called 

* 

4 

the  "Cadmium  Difference  Method". 

The  Trice  Method 

Threshold  and  fast  fission  type  foils  are  used  to  determine  the 
flux-spectra  in  the  fast  neutron  region  (Mev  region),  for  effective 
thresholds  of  presently  available  materials  range  from  approximately 
0.7  to  8  Mev. 

•  * 
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The  probability  of  finding  a  spectrum  with  the  general  character¬ 
istics  of  a  fission  spectrum  indicates  a  method  by  which  the  actual 
fast-neutron  cross  sections  can  be  transformed  into  idealized  thres¬ 
hold  types.  This  is  done  by  using  the  activation  formula,  Eq.  (2), 
to  determine  an  effective  threshold,  Eft. 

The  calculational  techniques  are  as  follows: 

1.  A  value  of  an  effective  cross  section  is  selected;  this  is 
denoted  as  <5r'  ,  a  constant.  (The  value  of  the  cross  sections  near  the 
effective  threshold  are  the  most  important  because  of  the  decrease  of 
flux  with  energy.) 

2.  The  activation  integral  is  equated  to  the  simplified  Integral* 


^°°  <T(E)  <j>(E)dE  =  (j>(E)dE 


-(00 


3.  (E)  is  assumed  to  have  the  form  of  a  fission  spectrum*. 

The  following  representations  have  been  proposed: 

(J>  (E)  d,  e_E  sinhVTa  (Watt) (5)  [l3*  ■]  or, 

(j)  (E)  cC  e"2/’965  5inhV2.29E  (thiv-  of  "California*) 

l  (b)  jj;  fctz.-t.7o3 


<f(E> 


values,  the  activation  integral  is 


ti.  Using  the  tabulated 
numerically  integrated. 

$.  The  value  of  the  simplified  integral  is  found. 


roo  f  00 

<|)(E)dE  =  Vj  CT(E)  -|>(E)dE 


*  These  functions,  normalised  to  unit  area,  are  tabulated  in  Appendix 
C. 
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6.  Ee  is  determined  from  the  fj)(E)  tabulation.  In  this  manner, 
the  true  cross  section  is  replaced  by  a  step  function,  the  threshold 
energy  being  determined  by  the  expected  spectral  form,  the  fission 
spectrum. 

As  an  example,  a  calculation  is  made  of  the  P31  (n,  p)  Si31 
effective  threshold. 


Q-  (3)  <j>(E)dE  -  23. 65 


(using  California  spectra) 


G  -  80  sb 


from  <|(E)  tabulation,  Efl  -  2.5  Her. . 

Figure  h  shows  the  transformation  made  in  the  1  (n,  p)  Sc  33- 
cross  section  by  this  method. 

A  list  of  the  effective  threshold?  and  saturation  cross  sections 
for  several  materials  are  listed  in  Table  1. 

In  an  actual  experiment,  the  foils  are  not  usually  irradiated  to 
saturation  and  are  not  counted  immediately  after  irradiation.  The 
following  formula,  derived  in  the  standard  manner,  gives  the  approxi¬ 
mate  value  of  the  integral  flux  from  the  experimentally  determined 
activity! 


oo 

|(e.)  --  ]  m 


>T 


d£  - 


h0  r 


I-  e 


-Jt 


(7) 
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Vhere:  S  =  activity  at  T  time  after  irradiation  in  dia/sec*mg. 
S'  =  microscopic  cross  section  in  cm2, 

H0  -  atans/mg  of  pure  isotope  present. 

/i  =  decay  constant  of  transmutation  product  or  .693/T(/l> 
$CE£)  c  \  <£>(£) dE»  the  total  neutron  flux  above  the 
energy  Ea— ■ the  integral  flux. 

(f>{£)  s  differential  neutron  flux,  see  Eg.  (2). 

T  =  time  after  irradiation, 
t  =  irradiation  time. 

£  s  counter  efficiency 

By  using  several  different  foils,  it  is  possible  to  make  a  plot 
of  the  integral  flux  versus  energy.  If  the  curve  drawn  through  the 
experimentally  determined  points  does  not  approximats  a  fiasion 
spectrum,  the  validity  of  the  spectral  knowledge  gained  is  uncertain. 
Another  defect  is  that  the  true  cross  sections  are  not  most  advanta~ 


geously  used  when  they  are  transformed  to  a  step  function. 

An  investigation  has  been  made  as  to  the  degree  to  which  this 
method  will  predict  the  correct  spectrum  in  the  absence  of  a  fission 
spectrum.  The  numerical  results  of  this  test  are  shown  in  ?ig,  5  an<* 
6.  The  spectra  assumed  were  as  follows* 

1.  ^)(E)  =  .0715  per  *5eT  for  E  less  than  14  Mev. 

(j)  (E)  =  0  for  E  greater  than  14  Mev#  (8) 

2,  <^(E)  =  *38  per  Mev  for  E  greater  than  0  Mev  and  less 


than  1  Mev. 
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Mkl 

Table  of  (c,p>.  (n.«0  ,  and  (n.  f)  Threshold 
Reactions 


1 5 


Reaction 

E. 

(Her) 

(harae) 

Half  Ufa 
of  Produat 

p3l  (q,)  S131 

2.5 

.075(aJ 

2.62  h 

s32  (np)  p32 

2.9 

.300^ 

14-3  d 

AI27  (np)  Mofl27 

5-3 

.080^ 

9.8  m 

Ni-58  (up)  C058 

5-0 

1.23 

72  d 

Si28  (np)  Al26 

6.1 

.190(d* 

2.3  a 

M^,24  (np)  Na2& 

6.3 

O 

O 

• 

15.0  h 

Al27  (n<0  Ka2^ 

6.6 

,110(o) 

13.0  h 

Wp237  (nf) 

•75 

1.32 

O2^8  (nf) 

1.45 

.54  M 

>  See  Text 

Th232  (nf) 

1-Z5 

.15  « 

J 

(a) 

Ref. 

(b) 

Ref. 

(c) 

Ref. 

(d) 

Ref. 

(e) 

Ref. 

3,  page  103. 

3*  page  108,  and  Ref-.  7.  page  106# 

2.  page  29. 

9,  page  248.  EEL-323  value  at  14  Her  used  for  normalization.* 

12,  page  343.  EEL-3?3  value  used  for  normalization. 
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<J>  (E)  ■  o8/B  per  Her  for  E  greater  than  1  Mev  and  lee# 
than  lit  Mev. 

ih  (E)  “  0  for  E  greater  than  lit  Mev.  (*}) 

3.  6  (E)  -  .0102(Uic£)  per  Mer  for  2  less  than  I4  Y*r. 

4>  (E)  =  0  for  E  greater  than  lit  Yet*  (t°) 

Ij.  (f>  (E)  =  .U5l9‘E//*965  w&rifuSZ  per  Kev  fcr  all  E 
(California  fission  spectrum). 

In  view  of  the  errors  encountered  in  making  any  foil  measurement, 
it  is  apparent  that  the  Trice  Methods  give  reasonable  results  f>3j0 
even  when  the  actual  spectrum  differs  vastly  from  the  fission  spectrum*. 
By  making  a  correct  assumption  regarding  the  true  spectrum,  it  is,  of 
course,  possible  to  calculate  the  thresholds  in  such  a  manner  as  t-e 
get  correct  results;  however,  the  labor  involved  in  doing  this  would 
not  be  insignificant. 

In  a  thermal  reactor,  a  spectrum  having  a  shape  combining  a 

fission  and  a  _JL_  spectrum  is  anticipated.  Fcr  this  type  of  situation, 
E 

it  appears  that  the  Trice  Method  would  give  cuite  acceptable  results. 

As  an  illustration  of  an  actual  reactor  measurement,  Fig.  7  la  a 
plot  of  a  portion  of  the  results  obtained  at  th9  Materials  Testing 
Reactor  in  195h»  The  prediction  of  this  data  is  that  the  spectrum 
measured  is  very  close  to  a  pure  fission  spectrum. 

The  fission  rate  (equivalent  to  saturated  activity)  of  fast- 
fission  foils  can  be  found  by  extrapolating  from  the  fission  product 
decay  curve  ,  cr  by  using  a  fa3t-ueutron  fission  chambe-. 


Point 3  are  the 
calculated  results 


j  Curves  aro  tho  true 
integral  spectra 

:rr 

zr-ut. 


5  6 

E  (Kev) 


Pi  rare  5 


frcniction  of  Spectra  oy  tho  Trico 


Joints  ruro  the 
oidenl  a  tod  results 


Curves  nro  tho  time 
integral  spectra 


PI  (J uro  6 


diction  of  In  outre. 
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In  neutron  fluxes  above  10*  neuts/cm^.sec,  only  the  first  technique  is 

r  ““ 

presently  suitable  |^10j16J.  G.  S.  Hurst  has  found  that  the  fission  pro¬ 
duct  gamma  decay  rate  (greater  than  l  Mev)  is  essentially  the  same  fcr 
plutonium-239,  neptuniun-237,  and  uranium-238,  and  is  apparently  inde¬ 
pendent  of  neutron  energy.  The  details  cf  this  analysis  are  not  pre¬ 
sented  here,  but  can  be  found  in  Ref .  $  and  6. 

.The  Cadmium  Difference  Method 

The  Cadmium  Difference  Method  utilises  the  resonance  activation 
integral  of  resonance  type  foils  to  achieve  spectral  knowledge  in  the 
l/K  region.  The  predominate  reason  for  this  is  that  the  resonance 
integral  is  better  known  than  the  true  cross  section.  The  discussion 
is  restricted  to  materials  having  resonances  at  energies  greater  ihir. 
100  ev.  The  resonance  absorption  by  cadmium  and  the  shape  of  the 
cadmium  transmission  curve  are  not  important  above  this  energy. 


r  00 

Resonance  Integral  =  \  T  dE  ...  (l  l) 

It  is  assumed  that  .CIO  inch  cadmium  covers  are  used,  an  isotropic 
flux  distribution  exists,  and  the  material  is  present  in  such  small 
quantities  so  as  not  to  perturb  the  flux  in  the  resonance  region. 

For  a  pure  l/v  activation  In  a  <E/E  type  cf  flux,  the  thermal 
cross  section  (2200  ir /set)  is  two  times  the  activation  integral. 

This  is  shown  by  the  follcvirgi 

CT  =  _i-cross  section  -  -r*i-  where  I  is  a  constant. 

▼  It 


data  taken  2.8  in.  from  end  of 
K3-3  bean  hole 


MTS  power:  30  Kf 


Fission  Spectrua 


Figure  7 


Fast-fteutron  Flux  Above  Threshold  vs.  Threshold 
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<Th. 


k  /or 

\  K"/=’  -* 
*•/ 


'  /{7T<  , 

— : -  -  2.c> 

z/or 

(12) 


Tha  total  epi-cadmiun  activation  is  proportional  to 


0_4£> 


% 


(o) 


-  vpsonar.c?  activation 
1/v  activation  ^  ^ 

When  the  foil  is  irradiated  without,  cadaiua  (tare)s  the  tctal 
activation  is  given  by  the  following* 

'00 

Abare  S  $th  rth4  ^Szl  (•*) 

Where s  3  the  satiated  activity  without  cadadux. 

/p  th  “  +-^e  tfc*fsal  r.eu^rcr.  flux  the  total  ther- 
nal  neutron  der.3it7  tines  2.2  x  10^  cis/st>c» 

-  the  22 OC  is/sec  activation  cross  section, 

y  o  “  4  constant. 

-  the  activation  :rcs3  section  aj  a  function 
of  energy. 

E^  -  the  erv-rgy  at  which  the  resonance  flux  j  elite 
tha  th'rzal  flux. 


^  "y/s 


The  lower  Unit  of  the  ab:ve  Integral  is  take::  as  Z^,  for  in  an 
actual  activation  sons  neutrons  are  absorbed  in  slowing  dcvn  from  ,k 
ev  to  thermal.  The  energy t  e\  for  which  the  dF./'Z  and  the  Maxwell 
distributions  are  equal  in  a  graphite  mode- rated,  room  temperature 
reactor  is  approximately  .2a  ev.  The  following  calculation  indicates 
the  importance  of  this  region  for  a  pure  l/v  foil. 


•:l  ^r< 


(II,) 


The  total  activation  cf  the  foil  irradiated  with  a  .010  inch 
cadmium  cover  in  an  isotropic  1/3  flux  is  as  feilevs* 


Acd 


JE 

T~ 


(17) 


Irradiated  without  cadmium*  the  total  activation  is  as  fcll-'v?* 

n,-4  ^  , 

A  =  1 4*-,  ,  i  f  ti?  *  4  \  <r4e 

'’Q  ^ty.  rk  T  <0  JfX 4  C  °  J  4  t 


The  cadmium  rati:;  C2,  is  gl7*r.  by  the  following* 


Qffi-9 


2k 


CR.  =  do/A,  i  -  I  +  i^±S±  *  ' ,z  ^Trk. 

o  Lt>  ,  rcs 

\  Sdsf£ 


CR-/  =  ^  ( ^h,  ♦  • 11 
•V 

Now,  (j>o  is  usually  much  less  than  6^  (  Ct  loo 

fore,  the  .12  <&0  ^  .oot  and  can  be  neglected. 


CZ-1  - 


St.  bat 


■•f 


Rearranging, 


**  4}t 


<j>  =,  -!*. 

0  (ct'i)  S"**'** 

w 


But, 


f  C-c/F/e  =  (ft«0  f  --/ir 


_  O*^) 


OO 


<2o). 

)}  there - 


( 2, ) 


(22) 
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The  usual  practice  in  making  a  measurement  in  the  1/i  region  Is 
to  determine  the  thermal  flux  by  cobalt,  gold,  or  indium  irradiation. 
The  cadmium  ratio  of  several  of  the  dete:tors  listed  in  Table  2  ia 
then  determined  by  irradiation  and  counting.  If  the  spe strum  ia  indeed 
the  l/Z  type,  the  determined  for  each  resonance  detector  should  be 
the  sauna.  The  differential  flux  ia  then  given  to  be  the  following* 


4>  m  -  &  !* 


(over  the  range  of  energy  determined 
by  the  detectors.)  (iV) 

i 

If  the  6afa  are  not  approximately  the  same,  little  can  be  said 
about  the  spectrum. 

In  all  such  measurements,  care  should  be  taken  to  insure  that 
the  correct  cadmium  thicknesses  are  used  (in  the  above  example  .010 
inches),  and  that  the  resonance  detectors  used  do  not  appreciably 
perturb  the  flux.  In  addition,  it  is  not  usually  advisable  to  irradi¬ 
ate  bare  and  cadmium  covered  foils  ia  the  same  location  at  the  same 
time. 

Discussion  and  Conclusions 


% 


By  combining  the  results  obtained  by  threshold  and  resonance  type 
foils,  a  crude  picture  of  the  neutron  flux-spectrum  of  a  reactor  can 
many  times  be  obtained.  A  1/Z  extrapolation  frcm  9100  ev  to  the  point 


at  which  it  neats  the  threshold  fell  curve  is  usually  mads  j  however* 
this  does  not  imply  measured  results  in  this  region.  At  the  low  energy 
end  of  the  spectrum,  the  spectrum  is  assumed  to  be  1^5  down  to 

Table  2 

Table  of  Eescnar.ce  Parameters 


Reaction 

Kn,r)I+1 

Resonance 

Energy  (ev)  ^'t 

Half  Life 

of  Product 

Cobalt-5? 

120  2.0 

- — 

5.3  7 

Manganese-55 

260  1.0 

•  1.33 

2.6  h 

Copper-63 

570  1.3 

.933 

4 

12.8  h 

Sodium-23 

1710  0 

.238 

15-  b 

Chlorine-37 

1800  o 

.079 

38.  m 

Vanadium-50 

3000 

.13U 

3.7U  a 

Aluminum-27 

9100  .68 

.159 

2.3  a 

•j*s  calculated  from  data  presented  in  Reference  J3U 


Using  020  in.  cadmium 
(.5  ev  cutoff) 

Experimentally  determined  using 
.010"  cadmium. 

not  be  used  until  checked  in  known 


•s  obtained  from  Reference  1C 


A,,  2-27<&th._ 

7  (cft-i)  U  + 


(cR-iXi^jy 

The  oC.’s  listed  above  should 


l/S  spectrum. 


17.  The  Role  of  Orthogonality 


Introduction 

In  the  following  sections,  several  flux-spectra  analysis  Methods 
are  discussed  which  require  a  minimum  of  spectral  assumptions  and 
which  use  the  activation  cross  sections  of  foils  in  their  true  form. 
That  is,  no  approximations  regarding  step  functions  are  made.  Before 
progressing  with  a  discussion  of  these  methods  in  detail,  it  is  best 
to  discus 3  them  in  a  general  manner. 

The  general  procedure  has  been  to  express  the  activations  as  a 
linear  combination  of  terms  which  depend  upon  some  parameter  describing 
the  spectrum.  The  number  of  the  equations  is  equal  to  the  ntndber  of 
foils  used. 

A  general  requirement  in  the  solution  of  systems  of  equations  is 
that  the  equations  themselves  be  independent.  Since  any  analytical 
method  of  obtaining,  flux-spectra  information  from  foils  will  be  depen¬ 
dent  upon  the  foil  cross  sections,  it  is  obvious  that  these  functions 
must  not  be  too  similar. 

Two  functions  are  said  to  be  orthogonal  if 


r  oo 


X.  (e)  X,(E)<3  «  C 

'O  1  J 

Wher.e,  I^(E)  and  Xj(e)  are  the  functions. 


(zs') 
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As  a  definition,  X^(E)  is  said  to  be  "independent"  of  Ij(E)  if 
the'  above  integral  is  true* 

If  ^°°  Xi1(S)X;j1(E)dE  =  0,  Xi1(E)  and  Iy^E)  are  said  to  be 

"dependent".  The  functions  1^  have  been  normalized  bo  that} 

00 
0 

When  the  value  of  the  orthogonality  integral  is  between  0  and  1 , 

1 

the  functions  are  said  to  be  partially  "dependent",  the  d6grse  of 
"dependency"  hinging  on  the  magnitude.  (The  orthogonality  integral, 
as  used  here,  is  really  the  inner  product  of  two  normalized  vector 
r  functions  in  Hilbert  space  —  the  magnitude  being  the  cosine  of  the 

apparent  angle  between  the  functions.) 

Threshold  Type  Foils 

In  practical  cases,  the  cross  sections  of  the  threshold  and  fast 

fission  materials  are  not  orthogonal*  That  is, 

00  v 
<T(s)  r(E)dE  \  0 

As  a  simple  illustrative  example  of  what  can  occur  in  solving 
nearly  dependent  equations,  the  following  is  presented! 

Let  two  hypothetical  parameters,  X  and  X,  be  represented  by  the 
following  equations! 

10X  +  2T  =  2U  =  A 


t 


51  +  21  =  £0  ■*  B 
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The  solution  is*  Z  *  U 

1  =  8 

Now,  let  there  be  a  four  per  cent  error  in  A*  Then, 

101  +  21  =  23 
91  f-  21  =  20 
The  solution  is*  1  =  3 

T  =  -  3.5 

This  result  shows  that  an  error  cf  less  than  5%  in  A  will  change 
the  calculated  value  of  T  by  more  than  $0%, 

Let  it  be  assumed  for  a  moment  that  Z  and  T  are  spectral  parameters, 
and  A  and  B  are  foil  activations.  Then  the  coefficients  of  Z  and  T 
must  be  related  to  the  activation  cross  sections,  for  the  cross  sections 
are  the  connecting  "links"  between  flux  and  activation.  Now  an  exper¬ 
imental  activation  error  of  is  certainly  not  unreasonable}  however, 
this  errer  would  generate  an  error  cf  ever  $0%  in  one  of  the  calcu¬ 
lated  spectral  parameters.  The  reason  for  this  is  that  the  coeffi¬ 
cients  are  nearly  proportional;  the  system  cf  equations  is  said  to  be 
"poorly  conditioned". 

The  preceding  discussion  implies  that  difficulty  may  be  experienced 
in  accurately  determining  spectral  parameters  from  any  system  of 
equations  in  which  two  cr  more  cross  sections  are  quite  "dependent". 

As  a  useful  technique  in  reducing  this  problem,  use  is  made  of 
the  following  identity** 

*  See  Appendix  D  fox  derivation. 
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a  =  J^°4(z5  •te'dS  5  ^ 

■Where,  ^(E)  =  ^  (the  integral  flux). 

Since  threshold  detector  cress  sections  increase  rapidly  near 
their  thresholds,  jL|L  is  largest  in  the  threshold  region  and  quit* 
small  elsewhere  (see  Fig.  3).  Using  the  alternative  activation 
representation,  it  may  be  possible  to  solve  a  system  of  equations  for 
<jjT  with  less  error  even  though  the  cross  sections  are  *dependent*,  for 
the  derivatives  of  the  cross  sections  may  be  quite  ■independent*. 

Since  threshold  type  foils  axe  used  in  the  fission  spectrum  region, 
and  the  flux  (either  integral  or  differential)  decreases  rapidly  in 
this  region,  it  is  convenient  to  introduce  the  following  relationship! 


( 28 ) 


Where,  ?  p  s  "  the  integral  fission  spectrum  function. 

*5  z  the  ratio  of  the  true  integral  spectrum 
to  the  integral  fission  spectrum.  This 
term  is  fairly  constant  in  energy. 

^  :  the  true  integral  spectnau 
Now,  the  activation  equation  can  be  written  as  the  followingt 


ry?FS 

Jo  x  F*S*  dS 
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Since  la  fairly  constant  in  energy,  a  set  of  threshold  type 
foils  can  be  considered  as  "independent"  if  t 

ito 

$*'  iSs  £Tj  «£•  -  o  (30) 

„  -Lp'i-  oZ  d£ 


for  all  j  detectors. 

The  above  equation  defines  orthogonality  of  the  cross  section 

derivatives  with  respect  to  the  parameter  Jr  ^  . 

'  F.S. 

Resonance  Type  Foils 

As  previously  discussed,  the  cros3  sections  of  resonance  type 
foils  are  characterized  by  resonance  peaks  superimposed  on  a  l/v  dis¬ 
tribution.  Since  the  resonance  peaks  cf  the  detectors  are  generally 
narrow  and  separated  in  energy,  it  might  be  supposed  that  the  cross 
sections  are  orthogonal.  Ihfortur.ately,  this  is  true  only  if  the  1 /▼ 
portion  can  be  eliminated  from  the  overall  activation.  If  the  's, 
Eq.  (lit),  of  several  foils  are  much  greater  than  unity,  the  l/v  acti¬ 
vation  can  be  neglected;  if  the  <£* 3  are  les3  than  unity,  the  major 
contribution  to  the  activation  is  due  to  the  l/v  cross  section.  Since 
most  presently  used  foils  have  significant  l/v  cross  sections  and  small 

SC  's,  they  are  highly  "dependent".  Sir.c9  the  integrals  of  the  alter- 

.  -  * 

nating  cross  sectional  derivatives  are  nearly  zero,  they  are  not  suit¬ 
able  for  making  integral  spectrum  measurements;  however,  they  are 
suited  for  making  integral  spectrum  measurements  if  '/"is  large. 

One  practical  solution  to  this  K  problem  is  discussed  in  Appendix  E. 
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Since  at  this  writing,  there  is  a  lack  of  good  resonance  cross 
section  data,  the  usefulness  of  any  analytical  method  utilizing  true 
cross  sections  is  sosewhat  academic.  If  and  wh6n  such  cross  sections 
become  available,  the  resonance  type  foil  may  become  a  very  useful 
device  for  measuring  the  differential  spectrum  in  the  epithermal  energy 
region. 


I 


we- f 

i 

t 

J  V-.  TRa  Wslygona?  *fethod 

i 

i 

* 

,  Basic  Philosophy 

The  basic  postulate  of  this  method  is  that  true  spectrm  para** 
meters  can  be  represented  by  a  linear  interpolation  between  a  fixed 
number  of  empirically  determined  spectrum  Tallies.  Figure  8  is  a 
graphical  representation  in  energy  space  of  this  postulate.  In  thia 
case  the  ‘s  are  taken  as  the  quantities  to  be  determined  by  acti¬ 
vation  and  calculation. 

Since  the  true  spectrum  is  probably  a  smooth  curve,  a  certain 
amount  of  error  should  be  expected  from  such  a  representation  in  energy 

«■  space.  By  transforming  the  energy  variable  into  a  more  appropriate 

variable,  it  should  be  possible  to  greatly  reduce  this  error.  Some 
possible  transformations  are  discussed  later  in  this  section* 

With  five  unknown  tj)  's,  as  indicated  in  Fig.  8,  five  different 
foils  with  independent  cross  sections  would  be  required  for  their 
determination  (five  unknowns  —  five  independent  equations). 

Application  to  Threshold  Type  Foils 

General  Equations  In  Energy  Space.  As  discussed  previously, 
threshold  type  foils  are  most  suited  for  integral  spectral  measurements 
in  the  Hev  region  (plutonium-239  is  an  exception  to  this  rule).  The 
following  equations  describe  the  polygonal  method  in  the  simplest 
mathematical  sense;  that  is,  all  correction  factors,  counter  geometries. 
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nuclear  densities ,  etc.,  are  assumed  to  be  included  in  the  saturated 
activity  terms  (Aa).  For  reasons  of  convergence  and  convenience,  the  - 
spectrum  is  assumed  to  be  exponentially  decreasing  above  Eg.  The 
energy  value  Ey  corresponds  to  the  highest  energy  spectral  parameter 
(the  energy  corresponding  to  the  ^  value  in  Fig.  8). 

The  General  Activation  Equation  is  as  follows 


Where s  A^  a  saturated  activity  of  i^  foil* 

4 1  -  differential  flux  at  energy  T  . 

-  differential  flux  at  energy  By 
s  activation  cross  section  of  ith  foil. 

4V  =  differential  flux  at  energy  B  . 

N  •  number  of  different  foils  used. 

For  N  foils,  there  will  exist  X  equations  of  the  above  type.  By 
algebraic  manipulation,  these  equations  are  transformed  into  the 
following* 


f'. 


»  I 
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In  matrix  syxdboliSM* 
(A)  *  (C)  (j) 
Solving  for  p  s 

U)  =  (C)'1  (A) 


(3  V) 
(3S) 


Since  threshold  detectors  are  most  suitable  for  determining  inte- 

-j" 

gral  flux,  Q  (2),  the  inverse  (C)  matrix  should  be  multiplied  by  the 
integral  matrix,  (I). 

(J)  =  (I)  <$)  =  (I)  (C y1  (A)  <3*3 

For  this  particular  polygonal  method,  the  (I)  matrix  is  the 
followings 


Note  that  all  elements  to  the  left  of  the  diagonal  are  zero. 

The  construction  of  the  integration  matrix  is  based  upon  the  following! 


J  -t,  + 


>  -  •  /.  3<j  4/J  (30) 

» 

Wher<,J  ^  j-1  Ej  =  2j  “Bj-1 
The  last  term  Is  l,39<j>jj,  for  ^  e-.72(E-Ejj)  <je  a  1.39 

In  summary,  the  following  procedure  is  indicated  for  obtaining 
integral  spectrum  information  by  this  method*  * 

1.  Several  foil  materials  whose  cross  sectional  derivatives  are 
quite  independent  are  selected. 

2«  (C)  is  computed  by  "hand".  A  set  of  E  's  are  selected  for  the" 

J 

integral  limits;  each  is  selected  from  the  energy  region  in  which  the 
ith  cross  section  is  varying  between  its  minimum  and  saturated  value 0 
(The  Ej*a  selected  by  using  the  polynomial  method  described  in  the 
next  section  may  be  the  best.) 

3.  (C)*A  is  computed  by  an  electronic  computer. 
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4*  The  (I)  (C)"1  matrix  is  computed  for  the  Ej's  used  in  the  (C) 
calculation.  An  electronic  computer  is  m03t  conveniently  used  for  this. 

5.  Upon  completion  of  the  nuclear  measurement,  (I)  (C)“^  is 
multiplied  by  the  empirically  determined  A^‘s  (a  simple  ■hand*  calcu¬ 
lation). 

6.  The  result  of  step  5  yields  the  ^  ‘s. 

7.  If  the  (I)  (C)"^  matrix  i3  poorly  conditioned,  a  new  set  of 
£^*3  are  selected  and  the  calculation  is  repeated. 

The  Definition  of  Gay.  Happy,  and  Sad.  As  mentioned. previously* 
a  spectrum  having  a  fora  similar  to  the  fission  spectrum  i3  expected 
to  be  found  in  the  Mev  region.  Thi3  fission  spectrum  function  is  a 
rapidly  decreasing  function  in  energy;  hence,  the  integral  flux  above 
8  Mev  is  much  less  than  that  above  2.  Since  the  polygonal  method,  in 
energy  space,  is  based  on  a  linear  extrapolation  between  energy  values, 
the  error  involved  in  approximating  the  true  spectrum  may  be  signifi¬ 
cant.  The  error  involved  is  the  integration  error  which  results  from 
approximating  a  continuously  varying  spectrum  by  chords  (see  Fig.  8), 

It  is  believed  that  better  results  may  be  obtained  if  a  new  function 
is  generated  that  will  give  a  linear  formulation  with  regard  to  the 
fission  spectrum,  for  the  fission  spectrum  is  a  first  approximation 
to  a  reactor  generated  fast  neutron  spectrum.  If  correctly  detemined, 
this  new  function  would  serve  the  same  purpose  in  the  Mev  region  as 
does  lethargy  in  the  1/E  region.  Like  lethargy  flux  in  the  1/E  region, 
the  differential  flux  with  respect  to  the  new  function  would  be  a 
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constant  in  the  Mev  region.  (The  differential  flux  In  energy  would  ha 


sharply  decreasing  with  increasing  energy.) 

It  has  been  shown  that  the  following  is  a  fair  approx imat ion.  of 
the  uraniun-235  fission  spectrum 

<f>  (E)  S  cVe~  Q-.775S  (3g)[ls66£) 

where  C  is  a  constant. 

The  new  function,  called  *gay* ,  is  defined  so  that* 

4>  (g)  *  a  constant  for  a  pure  fission  spectrum,  where  g  =  gay. 
The  derivation  of  gay  is  as  followsi 


<f>(e)de  =  err  c  d£  =  epe#)  <Jg 


(j)(g)  -  C  ff  e  775 F  de /jg  ~  consent  -  c, 


^  =  CZVT  c'77re  de 


-,ns£ 

ft  e  de 


let  y  "  E 


c/y  -  .TfS 


•l7Zfa 
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=  c, 


r  -7  ne* 

1  <>«•' 


let  u2-y 


2  u  Ju.  *  «*y 


The  first  term  of  the  above  equation  can  be  evaluated  directly} 
the  integral  can  be  evaluated  by  using  tables  of  the  standard  error 
integral. 

The  constant,  C^,  i3  arranged  so  that  &rj,T  is  1.0  —  this  corres¬ 
ponds  to  an  infinite  energy.  By  numerical  evaluation,  has  been 
calculated  to  be  equal  to  I.I3.  Figure  9  is  a  plot  of  gay  versus 
energy. 

The  relationship  between  $  (g)  and  £  (Z)  is,  by  definition, 
the  following  1 

$(E)  =  4(g)  -||-  -  '6(z)  cVe"  e“»775® 

=  ‘i£/(g)  *  Fission  Spectrua  C^O 
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Now,  at  any  gay,  the  differential  flux  can  be  found  by  multiplying 
0 ft  In  the  fission  spectra  tabulation  (Appendix  C)  by  the  appropri¬ 
ate  value  of  *  (g). 

For  reasons  of  presentation  clarity,  (g)  is  called  the  "sad* 
flnx>  (j>  (g)  dg  is  called  the  *  happy*  flux.  A  reason  for  these 

names  is  that  the  integral  of  a  poorly  1  behaved*  function  is  better 
■behaved*  |  hence,  'happy*  is  better  "behaved*  than  is  *aad*. 

General  Equations  in  Gay  Space.  With  the  exception  of  the  last 
tern,  the  polygonal  equations  in  gay  space  are  identical  in  form  to 
those  in  energy  space.  The  A^'s,  which  must  have  the  same  numerical 
value  as  before,  are  as  follows i 


Ai 


*  -6)1  qtp  dg 


ij  f  w  # 


(H 2 ) 


Note  that  the  Nth  term  is  of  different  form. 

Except  for  the  coefficients,  the  general  coefficients  are 
identical  in  form  to  those  used  in  energy  space. 


The  Nth  coefficient  is  as  follows* 


iLUiiiiiluu; 
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..  where  Us  *f.77l 
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The  general  procedure  for  utilizing  this  method  is  the  same  as 
that  used  in  energy  space.  Once  the  happy  flux  is  computed,  the  energy 
flux  can  be  found  by  multiplying  happy  by  the  appropriate  value  in  the 
fission  spectrum  tabulation. 

Discussion  and  Conclusions.  The  determination  of  the  proper  Ej's 
or  gj*s  is  probably  the  greatest  probJem  encountered  in  using  the  poly¬ 
gonal  method.  Calculations  using  hypothetical  cross  sections  have  been 
made  in  the  attempt  to  find  some  general  selection  criterion  whinh 
would  not  involve  numerical  iteration  procedures.  The  results  of  these 
calculations  have  not  been  successful;  however,  the  energy  parameters 
(Ej's)  found  to  be  suitable  for  the  method  described  in  the  next  section 
(the  Balynomial  Method)  may  be  the  correct  ones  to  use.  It  is  hoped 
that  this  can  be  investigated  in  the  near  future. 

If  cross  sections  existed  which  were  cf  the  idealized  threshold 
type,  then  the  proper  Ej's  or  gj/o  would  be  at  the  respective  thres¬ 
holds.  On  the  other  hand,  in  this  pleasant  situation  the  Trice  Method 
would  give  equally  accurate  results. 
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Application  to  Resonance  Type  Foils 
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V/hereas  gay  was  believed  to  have  general  usefulness  in  the  fission 
specvra  region,  it  is  believed  that  the  best  general  function  to  be 
used  in  1/E  region  is  lethargy,  , 

-  ln|a- 

where  E0  is  some  arbitrarily  selected  high  energy. 

In  a  puro  1/E  flux,  (tf.)  is  a  constant  and  is  related  to  d  ^»t) 
by  the  following t 

(It)  =  E?}(E) 

The  genefao.  equations  and  coefficients  of  the  polygpnal  method 
as  applied  to  the  1/E  region  are  of  the  same  form  as  those  used  in  the 
fission  spectrum  region,  When  used  in  the  1/E  region,  lethhrgy  replaces 
gay,  and  the  top  limit  of  the  last  integral  becomes  zero  if  the  calcu¬ 
lation  is  taken  to  the  Eq  of  lethargy.  It  would  probably  be  convenient 
if  were  selected  as  1  Mev  and  the  resonance  region  assumed  to  end 
at  this  energy,  (In  most  cases,  very  little  error  would  result  in 
doing  this.) 

If  the  limits  are  taken  in  ascending  order  of  energy,  the  A^'s 
cr  all  the  integrals  should  be  multiplied  by  a  negative  sign.  The 
reason  for  this  i3  that  lethargy  decreases  as  energy  increases. 

Since  the  resonance  type  foiJ.3  are  most  applicable  for  the  deter¬ 
mination  of  differential  flux,  the  integral  matrix  would  not  be  used. 

The  correct  u^'s  should  be  more  easily  determined  tLan  in  the  threshold 
detector  situation,  for  the  resonance  cross  sections  are  quite  sharply 


peaked.  If  the  ^'s  are  large,  the  u^'s  should  be  the  lethargies 
corresponding  to  the  maximum  cros3  section  values;  if  the  (£**3  are 
email,  the  u^'s  will  be  somewhat  different.. 

The  polygonal  method  is  not  presently  suitable  for  use  In  the 
resonance  region  because  of  the  lack  of  data  regarding  the  true  cross 
sections  of  suitable  materials.  However,  it  is  hoped  that  a  test 
calculation  involving  hypothetical  cross  sections  can  be  made  in.  th» 
near  future.  The  results  of  this  calculation  should  give  an  indication 
of  the  validity  of  the  polygonal  method  when  applied  to  the  1/E  spectral 
region. 


71.  The  Polynomial  Method 


Baalg.  aUlospphjr. 

The  basic  postulate  of  the  polynomial  method  is  that  the  differ¬ 
ential  flux  can  be  represented  by  a  polynomial,  the  degree  of  which  is 
one  less  than  the  number  of  detectors  used.  It  is  well  known  that  any 
function  can  be  represented  by  an  infinite  series;  however,  when  this 
series  is  reduced  to  a  polynomial,  large  errors  can  occur  in  the  repre¬ 
sentation.  In  order  that  representation  errors  be  reduced,  the  infinite 
series  must  be  rapidly  converging  over  the  interval.  If  this  ia  the 
case,  the  first  few  terms  of  the  series  will  be  sufficient  for  fairly 
accurate  representation  of  the  true  function. 

AppliceUca  la  IfrreshpU  TyE&  Eojlg 

General  Equations  in  Energy  Space.  The  assumption  is  made  that 
the  following  expansion  is  possible** 

(f)(£)  =  (j>0(? )  [«.  +  +  G'X5'’  *■  *  *  *  a/vv  £  J  6/6) 

*  If  one  of  the  foils  (such  as  plutonium)  has  an  appreciable  cross 
section  in  the  region  of  1/E  spectral  dominance,  the  following 
equation  should  be  used  in  place  of  Eq.  (46); 

4(f)  r  &0/e  f  4a(e)  f  t*'*e*’  f  "*  * e" '*  J 

The  representation  of  Eq.  (ij6)  is  not  satisfactory  because  the 
fission  spectrum  decreases  at  low  energies  whereas  the  1/B  spectrum 
increases. 
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Where:  N  is  the  number  of  different  detectors  used. 

the  normalized  fission  spectrum. 

The  activities,  Ai's,  are  written  as  follows* 

p®  CO 

A'l  ~  j  (fie)  <r^r)  de  =  ^  <fo(e)[&c  +  <t,e  +  *•*]  fye)  Jr 

ft  © 

Expanding* 


cm') 


o°°  c°° 

-  ao\  $>fr)  ri(e)  4*  f  ^  \  cf0U)  F  <rL(e)  Jr  + 

Jo  J© 

r® 

j  $0^  ^  *1(0  de  +  " '  (W) 

In  order  that  a  digital  computer  can  easily  handle  the  required 
calculation,  the  following  change  in  the  equations  is  made* 

(1)  I  is  substituted  for  E/10. 

r*  00 

(2)  \  4o0 r)  ic-  i •>  f't  I  ft  zero  — -  cJh.  djyroKiMxtt'fcK. 

WWer 

Now,  A.^  becomes  the  following* 


a"'  r-M 

Al  -  a,  l  (fjx)  c^(x)  oV  +  a,  \  ^6r)  X  (f*Cc)  Jx  t 

0  M) 

where  <£U)  =  10  C f>c(E) 

*  .  • 

and  <p(X)  -  <r)c(x)[  a0*  a(X  ►  a^  •  a^1  J 
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► 


Tot  H  foils,  each  with  relatively  *  independent*  cross  sectional 
derivatives,  the  following  matrix  equation  can  be  generated! 

(A)  =  (C)  (a)  (*o) 

Where  (C)  is  the  matrix  of  the  calculated  integrals;  that  isi 


C n  '  \  <r;(xi  J/ 

Jo 

c'"* 

CiZ=  \  X  <f>J.x)  3^  Or)  J/ 

e 

r*1  -h 

Cjj  -  \  xJ"*  4o(x)  <rt( t)  Jx 

o 

C1*  n  I 

\  * 

Jo 


Solving  the  matrix  equation  for  (a)» 

(a)  =  (C)’1  (A)  (S’/) 

Determining  the  a^'s  by  this  method,  <j)  (Z)  •  and  therefore  <£>  (E), 
can  be  theoretically  found  for  any  E  or  Z, 

Recalling  that  threshold  type  detectors  are  best  suited  for  inte¬ 
gral  measurements,  the  following  technique  should  be  used: 

Let«  (x)ty  -  (I)  (a) 


(fl) 
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where  (I)  is  the  integral  matrix  depending  upon  Z« 

Novi 

(I)  (a)  =  (I)  (C)-1  (A)  =  (  J  (X))  (*3) 

The  determination  of  which  energy  values  to  use  in  the  integration 
is  done  as  follows* 

1.  Several  X  values  are  selected  and  the  following  integration 
vector  is  computed  for  each  X.  (This  integration  vector  is  one  row  of 
the  integration  matrix*)  The  X  values  are  selected  from  the  effective 
threshold  energy  regions. 


( J)  1  (  [  <Pjx)J*  (  C  (ji  lt 1  '  J>t) 

''  *'■  1  (SI) 

* 

Upon  completion  of  the  above  computation,  each  integration  rector 
(J)  is  multiplied  by  (C)”1  to  form  (J)  (C)'1. 

For  one  particular  X,  the  resulting  vector  will  most  closely 
approach  one  of  the  following  forms* 


1. 

(  *. 

4  4  *  Uijj  «  •  4j 

z. 

(  «*« 

^5  . ) 

3. 

(  1 

tHjj 

Kj  . ) 

where 

the  hj  of 

each  vector  is  maximized  with  respect  to  the 

> 
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Form  1  is  approximated  when  X  values  corresponding  to  the  first 
(A1)  detector  are  used;  Form  2  is  approximated  when  X  values  corres- 
ponding  to  the  second  detector  (A^)  are  used,  etc. 

By  this  procedure,  one  determines  the  Ej  which  leads  to  the  best 
vector.  The  best  vector  is  the  vector  which  gives  the  minimum  error 
for  an  error  in  the  activation.  These  vectors  are  then  used  to  gener¬ 
ate  a  square  matrix,  the  ideal  form  of  such  a  matrix  being  the  following* 


(»>  -- 


:\ 


k.,  —  — 


(?S) 


Taking  the  product  of  this  matrix  and  the  matrix  gives  II  values 
for  (X).  In  this  manner,  it  is  believed  that  the  best  integral 
flux  information  is  generated. 


(j>  a))  =■  (2JU) 


For  the  ideal  case* 


<St>) 


$.  =  M.  = 
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Vhich  is  the  formulism  of  the  Trice  Method;  however,  is  this  .case 

the  true  cross  sections  have  been  vised,  and  the  assumptions  regarding 

the  spectrum  Lave  been  less  constraining* 

General  Equations  is  2az  Space.  The  polynomial  used  in  energy- 

space  was  selected  almost  arbitrarily;  however,  some  types  of  polynomials 

should  be  better  .than  others  when  applied  to  the  neutron  spectra  problem. 

A  theoretical  method  of  selecting  a  ■best*  polynomial  is  not  known  by 

the  author;  however,  a  polynomial  in  gay  space  seems  to  have  some  advan- 

« 

tages.  Let  such  a  polynomial  be  represented  by  the  followingi 

<f>W  o-o  *  *•  V, g""  <**> 

As  in  energy  space,  only  the  aQ  will  be  non-zero  if  the  real  spec¬ 
trum  is  a  fission ‘spectrum;  however,  in  gay  space  the  anticipated  spss- 
trum  is  more  properly  weighted,  for  the  flux -spectrum  assumes  more 
importance  at  low  energies  than  at  high  ones. 

A  practical  advantage  in  using  the  polynomial  method  in  gay  space 
is  that  the  calculations  axe  greatly  simplified.  The  loss  of  the 
^o(E)  weighting  factor,  particularly  in  the  integration  vector,  is 
responsible  for  this  simplification. 

The  procedure  for  using  gay  rather  than  energy  in  thd  polynomial 

i 

calculations  is  the  same  as  fcofore,  with  one  exception  —  S0(E)  is 


no  longer  used 
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An  example  of  this  change  follows  t 


LSI) 
{  60) 


The  integration  vector,  J,  becomes  the  followings 


(r)  =  (g  gW  g%  ■■■  8*/#)  (u) 

(This  can  bo  simply  calculated  by  hand  for  many  values  of  &•) 

'  Once  the  cross  sections  as  a  function  of  gay  are  computed,  the 
remainder  of  the  calculation  can  be  made  in  a  simple  manner,  In  cal¬ 
culating  the  cross  sections,  care  must  be  taken  with  materials  haying 
high-energy  thresholds,  for  the  gay  intervals,  a3  compared  with  those 
of  energy,  are  very,  small  in  the  high  energy  region. 

By  forming  a  matrix  of  the  best  (J)(C)"*  vectors,  the  happy  flux 
is  obtained,  as  before,  by  multiplication  with  the  activity  matrix. 

The  happy  flux  can  be  easily  transformed  to  energy  flux  by  simple 
conversion  of  gay  to  energy,  for  the  following  quality  is  true  by 
definitions 

J(g)=i£(E). 

where  g  and  E  correspond. 
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This  method  baa  not  yet  been  tested  using  materials  having  suffi¬ 
cient  independency  in  their  cross  sectional  derivatives.  It  has  been 
found  that  very  poor  results  will  be  obtained  if  two  of  the  cross  sec¬ 
tional  derivatives  are  quite  dependent}  i.e.,  sulfur-32  and  phosphor-31* 
(The  <^>(g)  matrix  becomes  very  poorly  conditioned.) 

Application  to  Resonance  Type  Detectors 

In  the  resonance  region  it  seems  plausible  that  a  polynomial  in 
lethargy  could  be  used.  A  possible  polynomial  is  the  following* 

(j>lu)  =  a0  +  a,  a  +  a2ux  *■  •••  u*~'  (63) 

Using  the  procedures  described  in  the  previous  subsection*  a 
matrix  equation  for  the  Ai's  can  be  found. 

(a)  a  (O'1  (A)  Litt) 

The  A^'s  calculated  in  the  foregoing  manner  are  the  coefficients 
of  the  lethargy  flux  polynomial?  hence,  in  theory  the  lethargy  flux 
becomes  known  for  all  lethargies.  Practically,  this  is  probably  not 
true,  for  the  polynomial  will  most  likely  represent  the  true  lethargy 
flux  at  N  lethargies  only. 

Each  of  the  N  lethargies  should  be  selected  by  the  same  method  as 
described  previously  for  the  threshold  type  foils.  Once  the  proper 
<£j'3  sure  determined,  the  approximate  lethargy  flux  throughout  the 
1/E  interval  can  be  obtained  by  constructing  a  smooth  curve  through 
the  (fj  points. 
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VII,  Numerical  Results  of  Test  Case 
Using  Polynomial  Method  in 
Energy  Space 

Is  an  indicator  of  usefulness  for  practical  work,  a  test  case  for 
the  polynomial  method  is  presented  in  this  section. 

The  test  is  made  using  the  following  spectra* 

1.  $  (E)  s  .0715  per  Mev  for  all  E  less  than  14  Mer* 

< p(E )  =  0  for  all  E  greater  than  14  Met.  (3) 

2.  (f>(E)  s  ,3 8/ET  per  Mev  for  E  greater  than  1  May  and  lass 

than  14  Mev. 

(p  (E)  -  ,38  per  Kev  for  E  greater  than  0  Mer  and  less 
than  1  Mot. 

(j)C E)  *  0  for  E  greater  than  14  Mer.  (f) 

3#  (p( E)  =  ,0102(14-2)  per  Mey  for  E  greater  than  0  Mer 

and  less  than  14  Her. 

(p( E)  =  0  for  E  greater  than  14  Mer,  Cto) 

4.  <£(E)  =  ,454e'E/*96-5  sinh  per  Mer  for  all  E, 

(The  California  fission  spectrum.) 

The  materials  and  specific  reactions  used  are  as  follovss 
Np237(n,f)  — C r 
U238  (n,f)  - - 
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<3 

* 

* 


s32(n,p)p32  — 

Si28(n.p)ll20  - 
Al27(n,d)Na2^  - 

These  materials  are  the  most  independent  of  the  group  of  seven  ' 
found  in  Appendix  B. 

The  (C)  matrix,,  as  calculated  .by  the  Wright  Air  Development  Center 
Univac-1103,  is  as  follows* 


78.6 

289.18 

103.15 

49.727 

-30.021  \ 

192.28 

88.39O 

32.836 

15.051 

-8.3735  | 

64.551 

26.782 

12.689 

£.8942 

1 

4.2864 

4.5335 

3.3009 

2.491 

1.9579 

1.6059  1 

.38225 

.35128 

3.314 

-.3211 

.31925/ 

The  inverse 

(C)  matrix. 

(C)"1,  is  as  follows* 

.009345 

-  .05221 

.09593 

-  .2460 

.446o\ 

.09481 

.6725 

-1.551 

4.090 

-8.467 

.03085 

-2.464 

6.658 

-20.11 

47.37 

.3889 

3-315 

-9.802 

34.85 

-94.07  / 

.1640 

-1-453 

4.538 

-I8.38 

57.35  / 
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(j)  (c)*1  vectors  for  several  energies  are  as  follows! 


E 

JCi 

JC 2  JC 

K 

*>•& 

7.569x10"** 

-6.979xlO“5  -3. 8 68x10-** 

-2.561x10-3  -3.398x10-3 

.7 

6.001x10“** 

5.755x10“**  -1.127xlD“3 

-5.552x10“**  -5.731*10-3 

.8 

1*. 633x10"** 

1.089x10“3  -1.599x10“3 

8.268x10“**  -6.955x10-3 

1.2 

1.011x10“** 

2.G25xlO'3  -l.UUxlO-3 

1.822xl0~3  ~U.582x10“3 

*1.!* 

7.999x10“6 

2.008x10-3  -6.072x10-** 

7.821x10-**  -1.152x10-3 

1.6 

4i.313xl0“5 

1. 800xl0“3  l*.  275x1.0-** 

-5.657x10-** 

1.621x10-3 

*2.8 

-2.282x10-6 

-l.U!i3xlO“**  U.571xlO“3 

-2.732x10-3 

2.031x10-3 

3.0 

1.327x10-5 

-3.259x10"**  U.  605x10-3 

-1.836x10-3  -2.753x10°** 

5.0 

2.177xl0‘5 

-2.1*31*10"**  i;092xl0“3 

6.7li6xlO"3  -l.077xlO“2 

5.2 

1.669x10-5 

-1.828x10-**  7.899x10-** 

6.821x10-3  -9.690x10^3 

*5.U 

1.193x10-5 

-1.278x10“**  5-302x10-** 

6.752x10“3  -8.352x10-3 

8.0 

-3*898x10-6 

3.923x10-5  -1.1*72x10“** 

9.351*10-1* 

9.1*93x10-3 

8.2 

-3.033x10-6 

3.015x10-5  -1.108xL0-** 

5.685x10-** 

9.983x10*3 

*8.1i 

-2.131x10-6 

2.082x10-5  -7.1ilhxl0-5 

2.101x10-** 

1.033x10-2 

Tho  starred  rotation  Indicates  the  best  vectors  to  use.  The 
selection  method  used  is  described  on  page  50*  It  should  be  noted 
that  certain  dependencies  and  poor  selections  of  energy  values  exist. 
The  best  vectors  are  summarized  by  the  following* 
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•6  Her  *  Ej  <■  .65  Hev. 

l.k  »£T^2<L5  Mev. 

2.1  Mev  <  Ej  3*0  Kev,  dependence  with  U-23S  and  SI- 2 3  la 
indicated* 

5.U  -'fev  <  E^,  lack  of  evidence  regarding  dependence  (should  find 
dependence  with  S-32)* 

8.U  Hev  C  E^,  again  with  lack  of  evidence  regarding  dependence* 
The  notation  us6d  above  indicates  the  energy  range  in  which  the 
beat  Ej  should  be  found.  The  vectors  can  be  optimised  by  performing 
further  calculations  within  the  energy  intervals.  It  is  noted  that  the 
Ej^  and  E^  guesses  were  quite  in  error. 

As  an  aid  in  visualising  the  qualitative  differences  between  the 
vectors,  the  vectors  can  be  multiplied  by  the  first  row  cf  the  (c) 
matrix  (this  generates  the  flux  equation  for  a  normalised  fission  spec» 
trum).  The  following  vector  array  has  been  computed  in  this  manner* 


E  (Mev) 

(J)  (C)-1  (A) 

*  .6 

R 

CO 

• 

".02 

-.025 

-.0116 

-.00130 

.7  . 

.705 

.163 

r.  072 

-.0025 

-.00219 

1.2 

.119 

.591 

-.093U 

.00825 

-.00175 

*  l.U 

.009U2 

.565 

-.0388 

•0035U 

-.000555 

1.6 

-.0508 

.526 

.0277 

-.00256  • 

.000616 

GNE-0 

E  (*»«*> 

(J>  (C)-1  <*) 

« 

55 

*  2.8 

-.00268 

-.01*21 

.296 

-.01235 

.000776 

3.0 

.0156 

-.095 

.297 

-.00830 

-.000105 

5.0 

.0256 

-.071 

.0701* 

.0301* 

-.001*12 

5.2 

.0196 

-.0531 

.0510 

.0308 

-.0037 

*  5.1» 

.OUjl 

-.0371 

CM 

.=» 

<«\ 

O 

. 

.0306 

-.00319 

8.0 

-.001*57 

.0111*5 

-.0095 

.001*21* 

.0362 

8.2 

-.00356 

.0088 

-.00715 

.00256 

.00381 

*  8.1* 

-.00251 

.00606 

-.001*79 

.00108 

.00391* 

From  the  preceding  array,  the  best  E  's  can  be  selected  (the 

3  # 

starred  values). 

By  taking  the  sun  of  the  elements  in  each  row,  the  fission  spec¬ 
trum  integral  flux  for  the  particular  Ej's  can  be  found.  The  results 
of  tills  addition  are  as  follows  < 


E 

<^(E)  calculated 

<^(E)  true 

.6 

.832 

.835 

•7 

.7  96 

.799 

1.2 

A 

l.U 

1.6 


623 

559 

5oo 


.625 

.560 

.501 
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g  (Me*)  3>(B)  calculated  ?(B)  true 


2.8 

•2l|0 

Os 

CM 

• 

3.0 

.209 

.209 

5.0 

.0513 

.0515 

5.2 

.01*1*6 

•01*14*. 

5.U 

.0386 

•0383 

8  JO 

.00521* 

.00521 

8*2 

.001*1*6 

.0014*5 

8.U 

.00378 

.00379 

It  Is  believed  that  the  results  can  be  considered  as  satisfactory* 

# 

Figures  10  and  11  are  the  graphical  results  of  the  test  calculation 
mentioned  earlier  in  this  section.  The  activities  used  were  numerically 
calculated  by  using  the  test  spectrum  and  the  true  cross  sections.  The 
(J)  (C)-1  vectors  used  were  these  corresponding  to  K  's  of  0.6,  l.b, 

2.8,  5«h,  and  8.U  Kev.  Worse  results  were  expected  from  the  5«U  and 
8.U  Mot  calculations  than  what  actually  occurred. 
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VIII.  Conclusions 

Several  methods  for  obtaining  neutron  flux-spectra  information 
have  been  presented.  Tne  simplest  of  these  is  the  Trice  Method; 
however,  this  method  is  also  the  most  inaccurate.  Where  the  Trice 
fast-neutron  method  gives  results  accurate  to  t  t-%  for  various  spec¬ 
tra,  the  polynomial  method  in  energy  space  reduces  this  error  to  i  10JC 
even  when  the  polynomial  integral  matrix  is  not  optimized.  (See  the 
Trice  method  and  polynomial  test  case  result’.)  The  polynomial  iwthcd, 
optimized  in  gay  space,  should  give  even  better  results.  On  the  other 
hand,  if  the  available  cress  sestiens  are  not  suitably  orthogonal, 
the  Trice  methods  are  the  most  useable. 

» 

|  The  most  difficult  method  to  work  with  is  the  polygonal  method, 

* 

for  Ej  selections  must  be  made  before  the  (C)  matrix  can  be  calculated. 

It  is  hoped  that  the  results  of  the  'polynomial  optimisation '  can 
he  published  in  the  near  future. 
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Appendix  A 


The  Calculation  of  Eirst-Collision  Energy  Absorption  from  Spectral  Data 


The  first-collision  dose  rate  in  a  material  undergoing  irradiation 


can  be  approximately  calculated  from  the  following* 

,  Co 


c  00 

D(fj)  -  '~r  \  e*raa>  4>(0  de  £ 


Me* 
See. 


Where*  ^  »  (■ ,  A  being  the  mass  number  of  the  material, 
A! C^Ce)  =  macroscopic  scattering  cross  section  of  the 
irradiated  material. 

E  a  neutron  energy  in  Mev. 

<$•(£)*  differential  neutron  flux  in  neutrons/ cm2.  sec-Mer 
0(Bj  )  =  dose  contributed  by  neutrons  of  energy  Bj  and 
greater. 

Knowing  <f>[ E)t  the  dose  rate  can  he  calculated.  Under  sane 
circumstances  only  values  of  dS  are  available.  If  the 

scattering  cross  section  is  approximately  constant,  the  dose  can  be 


calculated  ts  follows* 
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f:  J 


we, 


Integrating  by  parts  with  U~  £  ,  du=Je  t  d\f-  ^(e)d£ 
V=  §£<f>(£‘)j£'  : 

o 


P(fr- )  *.  r  jjV  {  £  -  Jto]]  [  -  ^  L§0  -}  (•>]  Jt  j 


•where: 

1  CO 


&  =  §  dE 


§Ctr)  r  J  <t>le')d£' 

§  4>(e')je'  -  J  4ce'i  te*  -  J  <j>le')de'  =  $o  -  ffe) 


After  substitution  of  limits,  the  following  is  obtained: 

0(e)  =  t  ^  f f. )  *  5*?f« JeJ 

If  £(E)  is  Jcrovn,  the  dose  received  in  srj  energy  interval  can  be 

obtained. 

For  an  energy  dependent  scattering  cross  section,  the  following 


formula  can  be  used 
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fj-  ?«j) 


vhera  it  hss  teen  assumed  thet 


-O 


dhe  derivation  of  the  above  formula  vss  made  in  tbs 


same  manner  as  fiat 


for  the  constant  scattering  cross  section, 
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Appendix  B 

Carnation  of  (n,p),  (n,4.),  and  (r,f)  Cross  Sections 

The  following  is  a  tabulation  of  cross  sections  for*  several  thresh  x.1 
and  fast-fission  type  detectors.  In  ussy  cases  the  cross-sections  h«7t 
been  "smoothed"  for  reasons  of  calculations!  convenience.  The  aorrcaa  of 
data  are  the  sane  as  those  listed  in  Table  1. 

Eev 


<n  -  Kp^7{n,f) 

<TZ -  t&8'n,f) 

<T3  -  p31(n.p)Si31 

<r4  -  S32(r.p)p32 

(Ts -  Al27(n,p)Mg27 

01  -  si23'-,?}^28 

0*7  -  Al27(E,d)5a24 

<q  -  Units  of  10“8  harna(hb) 

£  -  Ucita  of  Msv 
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E 

<r, 

<2 

«v 

<r„ 

<r 

0.1 

0 

0 

0 

0 

0 

0.2 

25 

0.3 

66 

0.4 

1^5 

0.5 

390 

• 

0.6 

640 

0 

0.7 

900 

1.5 

0.8 

1060 

s 

o.9 

1180 

11 

1.0 

1250 

15 

• 

1.2 

1350 

46 

1.4 

1370 

150 

1.6 

l4oo 

4l0 

0 

0 

1.8 

l4oo 

490 

5 

15 

2.0 

1420 

530 

10 

25 

2.2 

1450 

540 

28 

56 

2.4 

1450 

550 

35 

80 

2.6 

1450 

550 

50 

.80 

2.8 

1450 

• 

60 

100 

0 

3.0 

1450 

54 

130 

1 

3.2 

l46o 

71 

165 

2 

3.4 

1470 

550 

68 

220 

4 
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I 


\ 


€ 

<r. 

a 

*3 

<Tf 

fi 

r7 

3.6 

1470 

550 

71 

225 

6 

0 

0 

3.8 

1480 

80 

215 

6 

4.0 

1490 

80 

230 

7 

4.2 

1490 

80 

260 

8 

4.4 

1490 

80 

265 

10 

4.6 

1490 

80 

265 

16 

4.8 

1460 

80 

270 

18 

0 

5.0- 

1450 

,79 

275 

20 

10 

5-2 

1450 

79 

280 

25 

14 

5.4 

1450 

79 

280 

27 

27 

5.6 

1500 

79 

285 

32 

33 

5.8 

1600 

79 

285 

36 

44 

6.0 

1700 

79 

290 

4o 

60 

6.2 

1800 

79 

295 

44 

77 

0 

6.4 

1900 

79 

295 

48 

no 

1 

6.6 

2000 

79 

300 

51 

140 

3 

6.8 

2150 

79 

300 

55 

150 

4 

7.v 

2250 

79 

305 

60 

170 

5 

7.2 

2300 

79 

310 

62 

180 

7 

7.4 

2400 

79 

312 

70 

190 

9 

7.6 

2450 

79 

315 

74 

190 

13 

7.8 

2500 

79 

315 

76 

200 

17 

8.0 

2525 

550 

78 

320 

78 

210 

21 

1 


i 

J 

f 


caffi-9  72 


f 

<n 

<Ti 

<r» 

CV 

ft 

r. 

8.2 

2550 

550 

78 

222 

79 

220 

26 

8.4 

2550 

324 

80 

•H*  — 

8.6 

2550 

326 

37 

8.8 

2560 

320 

51 

9.2 

2600 

332 

58 

9.4 

2600 

334 

67 

9.6 

336 

75 

9.8 

338 

86 

10.0 

3k) 

93 

11 

350 

no 

12 

360 

• 

no 

13 

365 

no 

14 

2600 

550 

78 

370 

80 

220 

no 

9 
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A£P.endix  £ 

The  Normalized  Vatt  and  California  fission  Spectra 

The  following  is  a  table  of  ordinates  and  integrals  for  two  commonly 

used  uranium  -235  fission  spectra  functions*  The  functions  have  been 

Cm  I  rf*  . 

normalized  such  that  ^  (e)  d£  *  l  flW  j  <pc(£)  J?  */  - 

•  4^^  “  »483  ®  E  sinh  fix,  =  The  Watt  Spectrum 

1  -E/#g65  , _ . 

(pc(E)  =■  .454  e  sinh  v2V29fi  =  The  University  of  California 

Spectrum  £  •  3 

where  s.  B  is  the  neutron  energy  in  Mev* 


The  Vatt  Spectrum 


EjCMav) 

t(E) 

£EJ  cj>wCE)dE 

.00 

L  <4<««  -  $/• 

oa 

. 202 

^  .0140 

.9860 

0*2 

N 

O 

CM 

* 

.0377 

.9623 

0.3 

*306 

•0663 

.9337 

0*4 

.330 

.0987 

.9013 

0.5 

.345 

.1321 

.8679 
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BjCMsv} 

<j>w(E) 

0*6 

.352 

.1669 

.8331 

0.7 
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topendlx  D 

The  Derivation  of  an  Activation  Identity 


It  is  the  purpose  of  this  appendix  to  derive  the  following 
relationships 


A  -  f  <T(£)  dE  z  C  ~  Jr 

"a 


Where*  f\  =  the  activation  function* 

<rff)=  the  activation  cross  section  a3  a  function 
of  energy — threshold  type  foil. 

<f>(E )  5  the  differential  neutron  flux. 

£  Z  the  neutron  energy 

pOO 

j  r  the  integral  flux 

Integrating  the  first  integral  by  parts. 


CD  ^ 

,4*  fitf)  [f(o>-  ^  [g(o)  -  §(£) ] 


6e 


dr 


o  ~o 

CD 


Where*  |Y<?3-  \  Jg 

Noting  that  ^(oo)^Ot  «W  rCo)  *  Oj 

.A*  § 


^  4>Ce)  <r(£)  de  =  ^  ^ 


df 


Therefore. 
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Appendix  J 

Activation  Shielding  vith  Boron 


One  practical  solution  to  the  d  problem  discussed  on  pa^-e 
may  be  to  cover  seach  foil  with  a  sufficient  quantity  of  boron.  Boron 
has  an  almost  pare  1/v  cross  section  for  activation. 

Using  a  boron  covering  results  in  a  reduction  of  the  number  of 
low  energy  neutrons  traversing  the  foil.  Such  a  situation  could  be 
responsible  for  improving  results,  for  the  apparent  d's  are  increased. 

=  activation  due  to  resonance  cross  section 
activation  due  to  1/v  cross  section 

The  activation  shielding  with  boron  can  be  calculated  from  the 
following* 

/■"  -  £  I flfZr  i 

A  -  j  e  pee)  4l e>  c*e 

_t  i-rva-i 

Where*  e  #  *  the  boron  transmission  factor. 

=  the  boron  macroscopic  absorption 
cross  section  at  2200  m/sec  neutron 
velocity, 

E  =  the  neutron  energy  in  ev. 

t  ~  the  effective  thickness  of  the  boron 
cover,  (approximately  2  times  the 
true  thickness  in  an  isotropic  flux). 
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Q~(e)  r  the  actuation  cross  section  of  the 
resonance  detector* 

<}j(c)  -  the  differential  neutron  flux* 

The  above  equation  can  be  written  as  follows t 


A  -  \  ICO  J  +  <rVc^  <j,(e )  Ae 


r  r,eo“ 

~  T(e°'i  3  ^(e)  Je  4  \  r(e)  r  ^ 


/»  r  teats/ 

4>o  r(e0)  \  <rYcf  Ae/e  +  (  He)  <r  4(e)  Ae 


Vheret  T(E)  -  the  transmission  factor  previously  defined* 

<T  i*  the  total  activation  cross  section,  as  a 

t 

function  of  energy,  of  the  resonance  detector* 
trYei-  the  resonance  cross  section,  as  a  function  of 
energy,  of  the  resonance  detector. 

3tf)  -  the  value  of  the  transmission  factor  at 
the  energy  of  the  resonance* 
cf>o  -  the  neutron  flux  per  unit  lethargy. 

If  the  second  integral  of  the  preceding  equation  is  negligible 
compared  with  the  first,  an  opportunity  for  determining  (f>0  exists.  The 
second  tern  is  made  suitably  small  by  proper  election  of  t.  Thirty 
thermal  neutron  relaxation  lengths  should  be  sufficient;  therefore. 
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"t  should  be  approximately  .03  inches*  To  insure  that  thermal  neutron 
protection  is  present,  a  small  thickness  of  cadmium  should  be  placed 
betveen  the  boron  and  the  detector. 
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